US   NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CALIFORNIA 


8854 


RONEY 


1959 


THESIS 
R695 

Letter  on  cover: 

A  PILE  OSCILLATOR  FOR  THE  AGN  201 
REACTOR 


JAMES  R.  RONEY 


A  PILE  OSCILLATOR  FOR  THE  AGN   2C1    REACTOR 


***** 


James   R.    Roney 


A  PILE  OSCILLATOR  FOR  THE  AGN  2C1  REACTOR 


by 

James  R.  Roney 

Lieutenant,  United  States  Navy 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 
MECHANICAL  ENGINEERING 


United  States  Naval  Postgraduate  School 
Monterey,  California 

19  5  9 


■ 


A  PILE  OSCILLATOR  FOR  THE  AGN  201  REACTOR 

by 

James  R.  Roney 

This  irork  is  accepted  as  fulfilling 
the  thesis  requirements  for  the  degree  of 
MASTER  OF  SCIENCE 
IN 
MECHANICAL  ENGINEERING 
from  the 
United  States  Naval  Postgraduate  School 


ABSTRACT 

Pile  oscillators  have  been  successfully  employed  to  detect 
small  concentrations  of  in  purities  in  substances  such  is  boron  in 
graphite.   In  several  instances  the  technique  has  proved  more 
sensitive  than  the  commonly  used  "danger  coefficient"  method.  The 
danger  coefficient  method  using  the  Aerojet  General  Nucleonics  201 
Reactor  has  been  demonstrated  to  be  as  sensitive  as  the  pile 
oscillator  method  on  other  reactors;  thus  a  pile  oscillator  for  the 
AGN  201  Reactor  offered  the  possibility  of  a  further  improvement. 

Because  of  low  cost  and  simplicity,  a  design  utilizing  one  of 
the  installed  ion  chambers  in  the  reactor  ind  the  linear  amplifier 
on  the  control  console  was  chosen  for  construction  of  the  pile 
oscillator.  Additional  instrumentation  employing  the  otitput  of  the 
linear  amplifier  gave  excellent  performance  in  reproducing  a  test 
signal;  however,  when  the  pile  oscillator  was  put  into  operation 
satisfactory  reproduction  of  the  output  for  a  given  sample  could  not 
be  obtained.  It  is  felt  that  the  error  is  due  to  reactor  noise. 

It  is  concluded  that  a  pile  oscillator  utilizing  the  installed 
ion  chamber  and  linear  amplifier  of  the  AGN  2C1  Reactor  operating  at 
its  present  power  level  is  not  satisfactory  because  of  reactor  noise. 
It  is  doubtful  that  any  substantial  improvement  could  be  obtained 
unless  the  operating  power  were  increased  above  the  presently 
licensed  value. 
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1.   Introduction. 

The  time  behavior  of  a  nuclear  reactor  i.ay  be  affected  to 
varying  degrees  by  the  introduction  of  a  material  into  the  reactor. 
This  behavior  can  be  related  to  the  nuclear  cross  sections   of  the 
material.  Methods  have  been  devised  to  utilize  the  time  behavior 
characteristics  of  the  reactor  to  measure,  in  particular,  neutron 
absorption  cross  sections. 

Certain  important  impurities  have  microscopic  absorption  cross 
sections  which  are  very  large.  Suall  percentages  in  materials  will 
increase  the  total  macroscopic  cross  section  significantly.  Boron  is 
a  common  exai  pie  and  its  presence  in  graphite  or  iron  must  often  be 
quantitatively  determined.  The  methods  by  which  absorption  cross 
section  measurements  are  made  may  readily  be  adapted  for  measuring 
the  concentration  of  impurities.  This  is  usually  done  by  a  comparison 
with  a  standard. 

One  technique  is  called  the  danger  coefficient  method  < 1  [  .  A 
sample  of  material  is  placed  into  the  reactor  in  a  fixed  position, 
usually  the  center  of  the  core.  The  time  behavior  of  the  reactor  is 
noted  and  through  appropriate  relationships  a  comparison  made  with  a 
standard  sample. 

Another  method,  first  used  by  Langsdorf  12  \   ,  is  the  pile 
oscillator.  A  mechanical  drive  oscillates  a  sample  within  the  reactor, 
A  neutron-sensitive  ion  chamber  is  used  to  record  the  overall  power 
of  the  reactor.  The  oscillating  sample  causes  an  alternating 


Elementary  definitions  of  the  terns  and  concepts  which  are  taken 
fro:  nuclear  physics  and  reactor  theory  are  given  in  Appendix  I. 

2 
Numbers  in  brackets  refer  to  the  Bibliography. 


fluctuation  in  the  overall  reactor  flux.  The  overall  flux  containing 
the  variational  component  is  sensed  by  the  ion  chamber.  This  signal 
current  is  treated  so  that  the  variational  component  is  rectified  and 
integrated  with  respect  to  time.  The  output  voltage  of  the  integrator 
is  calibrated  with  a  standard  sample  permitting  the  use  of  the  device 
for  measuring  either  cross  sections  or  concentrations  of  impurities. 

In  addition  to  the  variation  of  overall  flux  caused  by  the 
oscillating  absorber,  there  is  a  significant  local  variation.  By 
placing  the  ion  chamber  close  to  the  absorber  the  local  signal  can 
also  be  detected.  If  thermal  neutron  cross  sections  are  being 
measured,  the  best  location  for  the  oscillator  and  ion  chamber  is  in 
a  region  such  as  the  graphite  reflector  where  neutrons  are  predominantly 
thermal.  In  this  manner  the  overall  reactivity  changes  induced  by  the 
sample  are  insignificant  compared  to  the  local  disturbance.  This 
technique  was  described  by  Hoover,  et.  al.  l3>  .  The  sample  was 
oscillated  through  an  annular  ion  chamber  located  in  the  graphite 
reflector  of  the  X-10  reactor  at  Oak  Ridge.  The  output  of  the  ion 
cha&ber  was  amplified  by  a  special  feedback  amplifier  which  passed 
only  the  variational  component  of  the  overall  signal.  This  output 
was  fed  to  an  integrator  after  appropriate  rectification.  This  method 
has  been  successfully  used  at  Bvtelle  Memorial  Institute  J4  1  . 

Comparing  the  pile  oscillator  method  to  the  danger  coefficient 
method,  Langsdorf  reported  an  increase  in  sensitivity  of  a  factor 
of  2C  over  danger  coefficient  measurements  with  the  same  reactor. 
Using  the  Aerojet  General  Nucleonics  201  Reactor  (AGN  201)  for 
danger  coefficient  measurements,  Carew  and  How  j  5  J  made  measurements 
as  precisely  as  had  been  done  with  other  reactors  using  pile 


oscillators.  It  therefore  appeared  possible  that  a  pile  oscillator 
for  the  AGN  201  might  produce  even  more  precise  measurements. 

It  was  decided  to  construct  a  pile  oscillator  for  the  AGN  201 
reactor  and  investigate  its  characteristics.  Many  possible  designs 
had  to  be  considered.  In  order  to  oscillate  a  sample  through  the 
core  of  the  reactor,  the  glory  hole  would  have  to  be  used.  The 
physical  construction  of  the  reactor  would  preclude  placing  a 
detector  close  enough  to  detect  the  local  flux  variation.  If  the 
oscillator  were  located  in  a  region  where  a  close-in  detector  could 
be  placed,  the  access  ports  or  thermal  column  would  presumably  be 
used.  In  this  case  the  overall  reactor  flux  would  not  be  greatly 
affected  and  the  reactor  would  oe  acting  primarily  as  a  flux  field. 
All  the  previous  devices  of  the  type  that  oscillated  within  a  thermal 

column  or  reflector  were  used  in  reactors  having  a  thermal  flux 

9  -2   -1 

of  the  order  of  10   neutrons  cm  "sec   .  Since  the  thermal  flux  of 

6  -2-1 

the  AGN  201  is  of  the  order  of  10  neutrons  cm   sec   ,  this  method 

would  be  at  a  disadvantage. 

Thus  the  type  of  oscillator  which  measured  the  changes  in  the 

overall  reactor  flux  was  indicated.  From  considerations  of  simplicity 

and  reasonable  expense,  a  design  was  base   on  the  employment  of  one 

of  the  permanent  ion  chambers  of  the  reactor  and  its  associated 

amplifier  on  the  Control  console.  It  was  felt  that  the  following 

procedure  would  be  appropriates 

1.   To  construct  a  mechanical  drive  mechanism  which  would 

oscillate  a  sample  of  material  within  the  reactor  with 

exactly  reproducible  motion  over  each  cycle  and  with 

accurate  variable  speed  control. 


2.  To  devise  an  instrumentation  system  which  would  reproduce 
accurately  an  input  which  simulated  the  actual  signal 
resulting  from  the  oscillation  of  a  typical  sample. 

3.  To  perform  experiments  oscillating  samples  within  the 
reactor.  Using  the  instrumentation  system  devised,  to  perform 
measurements  to  determine  if  the  output  reading  of  the  system 

,  is  reproducible  and  to  determine  the  sensitivity  of  the 

overall  technique. 

2.   Reactor  kinetics  for  an  oscillating  absorber. 

The  complete  formulation  of  the  eifects  of  an  oscillpting 
absorber  in  a  reactor  requires  extensive  developi  ent  and  involved 
expressions.  A  rigorous  treatment  is  available  from  Weinberg  I   6  [  . 
However,  since  the  pile  oscillator  method  to  be  used  with  the  AGN  201 
reactor  measures  overall  changes  in  neutron  flux,  the  problem  is 
considerably  simplified.  Since  interest  is  therefore  confined  to  the 
ti  ewise  variation  of  flux  rather  than  the  spacewise  variation,  the 
case  may  be  considered  in  which  the  time-dependent  part  of  the 
solution  of  the  governing  equations  is  separable  from  the  geometry- 
dependent  parts.  This  is  achieved  if  the  reactor  is  at  or  near  the 
critical  condition. 

In  the  following  treatment  the  assumption  of  a  homogeneous 

reactor  is  made.  The  one-group  space-independent  kinetic  equations 

1 
are: 


See  Appendix  I  for  development  of  these  equations  and  definitions 
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In  order  to  study  the  effect  of  an  oscillating  absorber,  let  it 
produce  a  change  in  k  which  is  sinusoidal  in  form.  Writing  this  as 

k(t)  =Ako  ejcJt   (j2  =  -1) 
and  limiting  ^k  to  SL;all  values,  A  k  «  1,  the  small  amplitude 
variations  of  n(t)  and  C.(t)  are 
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neutron  density  in  the  critical  state.  Substituting  these  expressions 

and  corresponding  derivatives  into  equation  (l)  gives 
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If  terms  (3)  and  (3)  are  neglected  as  being  second  order,  and  term  (4) 
is  neglected  as  being  small  in  comparison  with  term  Qp  ,  the  equation 

reduces  to 
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Again  neglecting  second  order  ten.  a  gives 
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Solving  this  equation  for  C,  and  making  the  substitution  into 
equation  (4)  gives 
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Discarding  the    last   term  under  the   assumption  Ak     n    «     n,      and 
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Thus  it  has  been  shown,  subject  to  the  approximations  made, 
that  the  ratio  of  the  variational  component  of  the  neutron  level 
to  the  average  neutron  level  is  proportional  to  the  input  change  in 
reactivity  for  a  given  value  of  CO,   W(jcj)  is  called  the  reactor 
transfer  function. 


The  results  of  a  set  of  calculations  for  different  parameters 
X  ,  and  for  varying  CO  is  given  in  reference  j  7J  .  The  particular 
results  of  interest  and  applicable  to  the  AGN  201  reactor  are  shown 


in  Figure  1. 
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Amplitude  of  W(jcj)  as  a  Function  of  cj 
Figure  1 


3.   Description  of  pile  oscillator  for  the  AGN  2(1  reactor. 

(a)   General  arrangement. 

Figure  2  shows  a  block  diagram  of  the  pile  oscillator 
constructed  for  the  AGN  2C1  reactor. 
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Figure  2 
Block  Diagram  of  Pile  Oscillator  for  AGN  201 

The  mechanical  drive  (l)  moves  the  sample  of  material  (2) 
through  the  core  of  the  reactor  with  uniform  linear  motion.  The  iotl 
chamber  (3)  detects  the  overall  neutron  flux  level  including  the 
change  induced  by  the  sample  and  supplies  a  signal  current  to  the 
linear  amplifier  (4)  installed  in  the  reactor  control  console.  The 
output  of  the  linear  amplifier  goes  to  the  first  stn.ges  of  the 
metering  syster.:  (3)  where  the  signal  representing  the  response  to  the 
perturbation  of  the  oscillating  sai  pie  is  removed  from  the  overall 
signal.  Tie  positive  part  of  this  voltage  is  integrated  with 
respect  to  time  by  the  integrator  (6)  froi.  which  the  output  may  be 
read  directly  on  a  voltmeter  or  recorded  on  a  strip  chart  recorder  (7) 
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for  graphical  analysis.  The  output  voltage  of  the  integrator  is  the 
criterion  for  the  comparison  of  an  unknown  sample  with  a  known  sample, 
the  frequency  and  number  of  cycles  being  fixed. 

(b)   Mechanical  components 

In  the  design  of  the  mechanical  drive  it  was  desired  to  provide 
an  arrangement  which  would  be  reasonably  simple  to  manufacture  and 
which  would  be  adaptnble  to  changes  found  necessary  in  the  course  of  «• 
the  investigation.  It  was  also  desired  to  drive  the  sample  in  such  a 
manner  that  the  wave  fonr  of  the  induced  change  in  the  multiplication 
factor  would  be  sinusoidal.  Since  the  absorption  icportance  function 
of  the  reactor  (see  Appendix  II)  can  be  approximated  by  a  cosine 
souared  function,  a  linear  ii  ovement  of  the  sample  is  equivalent  to  a 
sinusoidal  sar pie  motion  through  a  linear  flux  gradient.  Thus  the 
mechanical  drive  was  designed  to  produce  linear  movement  of  the 
sample  through  the  reactor  core. 

The  drive  mechanism  is  shown  in  Figure  3.  A  constant  speed 
electric  motor  drives  a  hydraulic  speed  reducer  to  which  is  connected 
one  sprocket  for  a  roller  chain.  Another  sprocket  is  mounted  so  as  to 
stretch  a  heavy  roller  chain  (3/4  in.  pitch)  between  the  sprockets. 
The  connection  of  the  chain  is  i.ade  by  a  tool  steel  pin  with  a  one 
inch  overhang.  This  pin  engages  a  vertical  brass  guide  slot  on  an 
aluminum  block  whic  h  is  free  to  rove  torizontally  on  two  guide  rods. 
Thus  the  travel  of  the  pin  causes  the  block  to  follow  in  a  linear 
motion  except  at  the  end  points  where  the  direction  of  motion  is 
reversed.  The  sample  holder  is  attached  to  the  sliding  block. 

The  sample  holder  is  made  from  3/4  inch  o.d.  aluninui  tubing 
having  a  wall  thickness  of  l/32  inch.  The  tube  was  filled  with 
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Figure  3 
Mechanical  Drive 
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reactor  grn.de  graphite  which  was  machined  to  close  diametral 
tolerances.  A  portion  of  the  aluminum  tube  was  cut  away  at  the 
midpoint.  An  interruption  in  the  graphite  at  this  point  permitted 
the  introduction  of  a  test  sample.  T'ie  graphite  filler  was  employed 
to  rtinii  ize  flux  variations  caused  by  the  scattering  of  neutrons  by 
the  sample. 

(c)   Instrumentation. 

The  instrumentation  system  converts  the  variations  in  neutron 
flux  into  useful  information.  The  signal  is  initiated  in  the  ion 
char,  be  r  in  the  reactor  and  amplified  by  the  linear  amplifier  on  the 
reactor  control  console.  Because  of  the  low  power  level  of  the  reactor, 

the  ion  chamber  current  at  rated  reactor  power  of  one  hundred 

_9 
milliwatts  is  only  of  the  order  of  If    arrperes.  The  linear  amplifier 

is  used  therefore  to  increase  the  power  level  signal  to  usable  levels. 

It  is  from  this  point  that  the  additional  instrui;  entation  takes  its 

inf om <  tion.  In  speaking  of  the  "signal"  the  reference  will  be  made 

to  the  output  of  the  linear  amplifier. 

The  signal  voltage  consists  of  the  sum  of  three  parts:  a  D.C. 

voltage  representing  the  average  reactor  power  level,  an  oscillating 

component  representing  the  response  to  the  perturbation  of  the 

oscillating  sample,  and  noise.  Since  the  oscillating  co.  ponent  is  one 

or  more  orders  of  magnitude  lower  than  the  D.C.  component,  it  is 

necessary  to  separate  it  fro»  the  overall  signal  to  permit  measurement. 

The  frequency  of  oscillation  is  quite  low  making  removal  of  the 

variational  component  by  usual  amplifier  techniques  difficult. 

An  amplifier  of  the  design  used  with  the  Oak  Ridge  oscillator  <3) 

was  constructed.  This  performed  -well  with  low  freouencv  sinusoidal 
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inputs  but  was  not  sufficiently  stable  when  the  inp\$t  voltage 
contained  a  significant  an.ount  of  noise.  This  was  eventually 
abandoned . 

The  problem  was  satisfactorily  solved  by  employing  operational 
amplifiers  of  the  type  employed  in  analog  computer  applications.  A. 
systet.  was  developed  which  subtracted  the  D.C.  level  fron  the  overall 
signal  leaving  the  variational  voltage.  This  voltage  was  rectified 
and  integrated  with  respect  to  time.  The  integrator  voltage  after  a 
selected  nuuiber  of  oscillations  was  the  basic  value  from  the 
rceasu recent. 

In  order  to  describe  the  instrumentation  system  in  detail  it 
will  be  necessary  to  first  describe  two  eleuentary  analog  computer 
circuits. 

A.  high  gain  operational  amplifier  is  symbolized  by 
If  the  output  is  fed  back  through  a  resistor  R„  to  the  input  and  a 
voltage  e,  is  fed  to  the  input  through  a  resistor  R.  : 


the  output  voltage  will  be 
R. 


e   a  - 
o 


R 


lj 


el  *  "  G  el 


where  G  ■  ~ 
Rl 


If  additional  voltages  are  fed  into  the  input  in  like  manner,  the 


negative  sum 


e  =  - 
o 


will  be  formed. 


If  li  =  R-,  =  RQ   etc,  G  is  constant;  the  shortened  symbol    — 
1     2,  o 

will  be  used  for  this  case. 

If  a  capacitor  of  capacity  C  is  used  for  feedbac':  in  place  of 
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the   resistor,         e 


,  sA^A-L       \] 


e      ,    the    integral    e   =-  e,    dt 

0  «U 


will  be  formed.  Likewise,  additional  inputs  may  be  included 

Hi 


and 


.-"  t 


e   = 


where   G 


R.C 

1 


If  all  the  input  resistances  are  e< ual ,  the  shortened  symbol 
will  be  used. 

It  was  stated  previously  that  the  D.C.  component  of  the  signal 
voltage  was  to  be  found  in  order  to  subtract  it  from  the  overall 
voltage  leaving  only  the  variational  component.  This  requires  a 
definition  of  the  variational  voltage.  This  may  be  done  in  the 


following  manner: 


t 


var 


dt  6: 


var 


max 


That  is,  the  time  integral  of  the  variational  voltage  will  not  grow 
with  time  but  will  remain  within  the  limit  given.  This  may  be 
expressed  functionally* 


e 
-e 


<t>- 


-     I  f  (e  "  e  )  dt  ± 
RCJo 


var 


RC 


max 


where  e  is  the  total  signal,  and  e,  is  the  instant meous  value  of  the 
D.C.  component.  In  other  words,  if  the  integral  of  (e  -  e,  )  is  formed, 
the  value  of  the  integral  will  not  grow  with  time  but  will  remain 
within  the  limit  given. 

With  this  idea  in  mind  the  following  scheme  suggests  itself 


e   ~  -e. 
o      b 
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Although  the  behavior  of  this  system  can  be  described  in  precise 

terminology,  it  nay  be  understood  by  considering  the  effect  if  e   is 

o 

not  eoual  to  the  defined  e,  .  If  e   is  nui.  erically  too  small,  there 

bo  j  t 

will  be  a  net  positive  voltage  input.  This  will  generate  a  larger  e 
but,  since  this  is  negative  by  the  sign  reversing  characteristic  of 
the  integrator,  the  net  positive  input  will  be  diminished  until 
equilibrium  is  established. 

Thus  far  no  account  has  been  taken  of  the  possibility  that  the 
average  power  level  of  the  reactor  may  be  drifting.  This  would  have 
the  effect  of  making  e   drift.  This  can  be  minimized  by  selecting 
proper  values  of  resistance  and  capacitance  in  the  integrating 
an.plifier  so  that  it  will  follow  the  change  in  the  average  voltage. 
Generally  speaking,  the  shorter  the  time  constant  of  the  integrating 
ifier,  the  better  it  will  follow  the  drift.  However,  if  the  time 
constant  is  too  short,  it  will  tend  to  follow  the  oscillating 
component  of  Die  input  as  if  it  were  a  drifting  D.C.  Thus  a 
compromise  is  reached  based  on  the  maximum  drift  anticipated. 

With  the  D.C.  component  determined,  the  variational  component 
can  be  found:  

-a ^s e    =  e  -  e, 

var        b 


In  order  to  render  the  effect  of  drift  insignificant,  two  more 

amplifiers  was  added  repeating  the  above  loop.  The  variational 

voltage  is  then  rectified  by  a  diode  leaving  only  the  positive  halves 

which  is  the  information  signal,  e.  „  .  The  use  of  half  wave 

to      info 

rectification  insures  that  the  part  of  the  signal  which  becomes 

e.  „   corresponds  to  the  sample  position  being  in  that  part  of  the 
info 
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reactor  where  the  cosine  waveform  is  closely  approximated  (see 


Appendix  II).  The  rectified  voltage,  e.  .  ,  is  fed  as  an  input  to 


info 


an 


integrator, 


var 


£ 


u 


info 


e.  p  dt 

info 


I 


\os ^x 


-jT 


kn 


In  order  to  check  the  reproducability  of  this  system,  a 
simulated  input  was  formed.  This  consisted  of  an  electronic  function 
generator  with  a  battery  in  series.  This  gave  a  D.C.  voltage  with  a 
sinusoidal  variation  superimposed.  The  magnitudes  and  frequencies 
were  chosen  to  correspond  to  the  output  of  the  linear  amplifier.  The 
system  was  allowed  to  reach  equilibrium  such  that  it  balanced  the  D.C. 

corpcnent  and  then  integrations  were  performed  for  a  fixed  number  of 

I 
cycles.  This  process  was  repeated,  and  the  resulting  output  voltages 

of  the  integrator  were  read.  The  readings  were  always  the  s:>me  within 

a  fraction  of  a  percent.  A  slightly  drifting  voltage  was  added  with 

no  noticible  change. 

A  schematic  diagram  of  the  amplifier  system  is  given  in  Figure  4. 
Illustrations  of  the  physical  arrangement  are  shown  in  Figure  5. 

A  strip  chart  recorder  was  included  in  the  instrumentation  so 
that  the  final  integrator  output  or  any  intermediary  voltage  could  be 
recorded. 
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Figure  5 

Instrumentation  Equipment 
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4.   Operation  of  the  pile  oscillator. 

In  order  to  make  a  measurement ,  a  standard  specimen  of  iron  was 
inserted  in  the  sample  rod  and  the  rod  positioned  in  the  drive 
mechanism  so  that  the  sample  was  in  the  center  of  the  reactor  core 
when  the  drive  mechanisi  was  in  its  i  id-position.  The  amplifiers  in 
the  instrument  system  were  balanced  and  all  inputs  and  condensers 
shorted.  The  drive  mechanism  was  started  and  the  reactor  control  rods 
adjusted  until  the  trace  on  the  linear  recorder  on  the  control 
console  showed  no  tendency  to  drift.  This  required  several  minutes. 
A.11  shorts  were  removed  from  the  instrumentation  system  except  for  the 
final  integrating  capacitor.  After  several  minutes  the  instruments  had 
steadied  out  and  an  integration  could  be  commenced.  This  was  done  when 
the  drive  mechanism  was  in  an  extreme  "out"  position  and  continued  for 
a  fixed  number  of  cycles  after  which  the  output  of  the  integrator  was 
recorded.  The  integration  for  the  same  sample  was  repeated  for 
approximately  five  readings.  Then  the  sample  was  removed  and  replaced 
by  a  sample  having  an  unknown  percentage  of  impurity  and  the  procedure 
repeated.  In  this  manner  the  unknown  sample  could  be  related  to  the 
standard  sample  which  was  for  all  purposes  free  of  the  impurity. 

Several  difficulties  are  inherent  in  the  system  which  limit  the 
total  time  allowable  for  the  measurement  of  each  sample.  It  is 
difficult  to  maintain  the  average  level  of  the  reactor  at  precise 
criticality.  The  nature  of  the  instrumentation  system  is  such  that  it 
must  be  balanced  at  least  every  half  hour.  The  choice  of  component 
values  in  the  instrumentation  was  to  a  large  part  governed  by  the 
ability  to  compensate  for  drift  effects.  This  in  turn  required  that 
the  oscillating  frequency  should  not  be  too  low  or  the  instruj  ents 
would  Lend  to  follow  t  ,e  oscillation  as  if  it  were  D.C.  drift  rather 
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than  the  information  bearing  signal. 

The  choice  of  oscillating  frequency  must  take  other  factors 
into  account.  From  the  plot  in  Figure  1  it  would  appear  that  a  lower 
frequency  of  oscillation  is  desirable.  However,  statistical 
considerations  raise  the  possibility  of  advantage  from  higher 
frequencies  since  a  greater  number  of  observations  can  be  accumulated 
in  the  same  total  time.  This  argm,  ent  assues  that  the  relative  error 
for  each  cycle  does  not  increase  as  the  frequency  of  oscillation  is 
increased.  Weighing  all  the  contrasting  factors  it  appears  that  a 
frequency  of  oscillation  near  the  inflection  point  of  the  curve  in 
Figure  1  gives  the  best  compromise.  This  corresponds  to  a  period  of 
about  two  seconds  for  each  cycle.  It  should  be  noted  that  a  complete 
cycle  of  the  drive  mechanism  produces  two  cycles  of  reactivity 
oscillation.  A  set  of  measurements  was  made  on  the  same  sample  at 
varying  frenuencies  of  oscillation  but  for  the  same  total  times.  This 
confirmed  the  cycle  time  of  about  two  seconds  to  be  proper. 

In  order  to  test  the  accuracy  of  the  systei  ,  a  set  of  samples 
,  was  chosen  which  had  varying  amounts  of  boron  as  an  impurity  in  iron. 
These  samples  had  been  measured  for  impurity  both  by  spectrographs 
methods  and  by  the  danger  coefficient  method  in  the  course  of 
investigation  reported  by  Carew  and  Mow  l  5  >   .  A  typical  group  i 


is 

tabulated  in  Table  1. 


Table  1 
Samples  of  Iron  Containing  Boron 
Sample   Weight  (grams)   Percentage  of  Boron   Reactivity  U*f  ) 

200 
210 
232 
218 
267 
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A 

6.381 

0.0001 

B 

3.351 

0.0020 

C 

6.37C 

0.CC50 

D 

6.239 

0.C072 

E 

6.261 

0.0250 

5.   Discussion  of  results. 

(a)   Magnitude  of  error. 

In  order  to  discuss  the  error  of  pile  oscillator  measurements, 
a  sample  standard  deviation  is  employed  as  a  basis  for  determining 
the  sensitivity  of  the  technique.  The  sample  standard  deviation  will 
be  here  defined  as  the  standard  deviation  of  a  single  integration  as 
determined  from  several  consecutive  integrations  with  the  same  m 
material  sample.  The  sample  standard  deviation  can  be  projected  to 
determine  the  sensitivity  of  a  measurement  for  which  the  integration 
process  is  repeated  wit  i  the  same  material  sample.  The  standard 
deviation  of  the  arithmetic  mean  of  the  integrations  would  diminish 
by  the  square  root  of  the  number  of  repetitions.  A  set  of  typical 
readings  for  one  sample  is  shown  in  Table  2  demonstrating  the  idea. 

Table  2 

Integrator  Readings  from  Integrations  on  Sample  E 
Cycle  time:  4.5C  seconds;  50  cycles  each  integration 

Reading  (volts) 


29. 

.3 

26, 

,2 

28, 

,c 

29, 

,7 

25, 

,9 

Average  =27.8 

Sample  standard  deviation  =  1.7  or  6% 
Reactivity  of  sample        270yup  appro*. 
Magnitude  of  error  16a»^>  approx. 

The  error  of  a  measurer ent  would  be  diminished  by  5 
approximately  if  five  integrations  were  made  on  the  same  sample. 
The  result  of  numerous  measurements  such  as  the  one  tabulated  was 
that  the  sample  standard  deviation  for  the  method  corresponded  to 
an  error  of  from  20  up  to  10  up  .  In  order  to  reduce  this  by  a  factor 

20 


of  ter,  would  require  deten  ining  the  arithmetic  mean  from  about 
ICC  integrations  for  eacii  sample.  If  each  integration  is  as  rapid  as 
two  minutes,  over  three  hours  would  be  required  for  the  measurement 
on  a  single  sample.  External  influences  would  undoubtedly  come  into 
play  in' this  time.  Employing  the  danger  coefficient  technique  in  a 
repetitive  manner  over  a  corresponding  elapsed  time  with  the  AGN  201 
reactor  gave  a  sensitivity  of  1  /up  . 

Since  the  instrumentation  syster  was  capable  of  accurately 
reproducing  the  test  signal,  examination  of  the  experimental  input 
signal  for  variations  is  indie  ted.  The  oscillating  signal  which  is 
the  input  to  the  integrator  was  recorded  on  a  strip  chart  recorder, 
an  exact  tracing  from  a  typical  run  is  shown  in  Figure  6. 


A 


Figure  6 
Integrator  Input  Signal 

The  severe  variations  are  apparent.  These  variations  are  a 
consequence  of  the  total  noise  in  the  system.  A  discussion  of  this 
noise  is  desirable. 

())      System  noise. 

The  total  noise  appearing  at  the  final  integrator  input  will  be 
called  the  system  noise.  Noise  is  simply  defined  as  unwanted  signal!  8 >, 
The  system  noise  may  be  traced  to  three  principal  sources:  reactor 
noise,  ion  chamber  noise,  and  amplifier  noise. 
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Although  the  equations  of  reactor  kinetics  are  written  for 
continuous  functions,  the  actual  physical  processes  occur  in  discrete 
steps.  At  sufficiently  high  power  levels  the  number  of  events 
occuring  are  sufficient  to  limit  deviations  from  the  most  probable 
pattern  of  events.  When  low  powers  or  small  perturbations  are 
considered,  random  fluctuations  become  significant. 

The  reactor  noise  results  from  fluctuations  in  neutron  flux 
which  are  caused  by  the  discrete  nature  of  the  neutron  creation 
process.  It  was  noted  that  the  delayed  neutrons  come  from  radioactive 
decay  processes  and  are  therefore  subject  to  fluctuations.  Both  the 
concentration  of  delayed  neutron  precursors  and  their  decay  times  are 
subject  to  random  variations.  It  is  convenient  to  express  the  change 
in  flux  which  is  caused  by  the  oscillating  sample  in  the  fori  of  a 
comparison  with  the  variations  in  flux  not  caused  by  the  sample 
(i.e.  reactor  noise).  This  is  simply  a  signal  to  noise  ratio. 

Frisch  and  Littler  S9>  derived  an  expression  for  the  ratio  of 
signal  to  noise  for  a  single  observation  of  an  oscillating  sample: 


S(6_>)  =  (Ak) 


2YNf  T  A"(<o) 


where 


1  +TA  (w) 
(^k)   =  change  in  multiplication  produced  by  the  sample 


T   s  probability  of  the  detector  detecting  a  fission 

\ 


X, 


-   mean  fission  rate 
T      =  total  time  of  the  observation 
A.((jj)      =  amplitude  response  of  the  reactor. 
Rewriting  this  as    S  (cJ )  =  (^k)  ~/2Nf  T  G     where 


G 


A  (o) 


_1  +  A"(CJ)^ 
shows  the  signal  to  noise  ritio  with  respect  to  reactor  noise 
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improves  as  the  square  root  of  the  ratio  between  two  power  levels. 

Ion  chamber  noise  has  been  demonstrated  experimentally  -T 1 C J  to 
vary  inversely  as  the  square  root  of  the  ion  chamber  current.  Since 
current  is  proportional  to  neutron  flux,  this  varies  in  the  same 
manner  as  the  reactor  noise.  Tims  the  total  noise  appearing  at  the 
linear  amplifier  varies  inversely  as  the  square  root  of  the  reactor 
power  level. 

The  third  type  of  noise  is  that  produced  in  the  electronic 
equipment.  This  was  found  to  be  primarily  sixty  cycle  hum.  This  was 
attenuated  in  the  measuring  scheme  employed. 

The  precision  oT  a  measurement  using  this  technique  will  be 
es  entially  determined  by  the  system  noise  described.  Since  this 
increases  as  the  power  level  of  a  reactor  is  reduced,  it  will  be  a 
formidable  obstacle  to  the  success  of  a  pile  oscillator  on  a  low 
power  reactor  such  as  the  AGN  2C1. 

In  reference  {llj  an  expression  was  derived  which  is  similar  to 
that  of  Frisch  and  Littler  discussed  previously.  In  this  case, 
account  was  taken  of  the  correlation  involved  in  the  type  of 
experiment.  Tiie  expression  for  signal  to  noise  ratio  is 


—•  =  "  ~\j    — - —    with  F  the  mean  fission  rate  of  the 
reactor,  T  the  total  time  of  the  observation  on  the  given  sample, 
and  ^  a  factor  involving  the  mean  number  of  neutrons  released  per 
fission.  For  a  typical  measurement  (a  single  integration),  this 
expression  gives: 

F  =  f.l  watts  x  3  x  10    fissions/watt  sec. 

9         / 
=  3  x  10   fissions/sec. 

>  N 

T  ■  2.5  x  10   sec.  (a  typical  value) 


)5Ak 


|  -  105Ak 


Thus  for  a  signal  to  noise  ratio  of  unity,   k  must  be  of  the 
order  of  10  '  which  is  approximately  ten  roicro-rho.  It  would  appear 
that  the  limitation  on  the  results  has  been  placed  largely  by  the 
reactor  noise. 

The  effect  of  the  noise  on  the  danger  coefficient  method  is  slight, 
In  the  course  of  an  observation,  the  noise  is  involved  only  when  it  is 
necessary  to  determine  the  reactor  period.  If  the  criterion  for  this 
determination  is  the  position  at  which  the  trace  of  a  linear  amplifier 
output  on  a  strip  chart  recorder  crosses  a  given  ordinate,  the  average 
value  of  the  flux  level  can  be  read  through  the  noise  by  drawing  in  a 
line  which  is  the  best  estimate  of  the  average  signal  and  accurately 
observing  where  this  line  crosses  the  ordinate.  A  human  operator  has 
considerable  ability  in  finding  the  average  value  of  such  a  signal 
through  noise.  This  technique  is  illustrated  in  Figure  7. 


Intersection  accurately  determined 


time 


Figure  7 
Time  Determination  from  Strip  Cha,rt  Recorder  Trace 

Thus  the  statistical  advantage  of  a  pile  oscillator  over  danger 
coefficient  observations  is  open  to  considerable  doubt.  The  noise 
factor  may  be  so  disturbing  as  to  outweigh  all  other  considerations. 
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6.  Conclusions. 

It  is  concluded  that  a  pile  oscillator  for  the  AGN  2C1  reactor 
utilizing  the  installed  ion  chamber  and  linear  amplifier  does  not 
give  satisfactory  results.  It  is  not  possible  to  attain  with  the 
equipment  employed  the  precision  obtainable  with  danger  coefficient 
measurements  on  the  same  reactor  (and  using  less  eouipi  ent). 

Since  the  instrumentation  system  consistently  reproduced  a  test 
input,  output  variations  with  the  experimental  input  are  attributed 
to  the  nature  of  the  information  contained  in  the  input  signal.  The 
reactor  noise  appears  to  be  the  major  cause  of  the  variations. 

In  order  to  improve  the  sensitivity  of  the  pile  oscillator, 
either  one  or  both  of  two  things  must  be  done.  The  quality  of  the 
input  signal  could  be  improved  by  the  use  of  a  better  ion  chamber 
placed  closer  to  the  core.  Placement  close  to  the  path  of  oscillation 
might  permit  detection  of  the  local  signal  to  increase  signal  to 
noise  r\tio.  Another  improvement  would  result  from  operating  the 
reactor  at  a  higher  power  level  than  the  presently  licensed  limit. 
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APPENDIX  I 

ELEMENTARY  BACKGROUND 

In  a  nuclear  reactor,  the  reaction  sequence  by  which  neutrons 
are  produced  and  lost  is  self  sustaining  if  proper  conditions  exist. 
This  is  described  by  the  multiplication  factor,  k,  being  unity.  This 
is  the  "critical"  condition  for  the  reactor.  If  k  is  greater  than 
unity  the  reactor  is  supercritical  and  the  neutron  population  will 
increase.  Conversely,  if  k  is  less  than  unity,  the  reactor  is 
subcritical  and  the  neutron  population  will  decrease. 

Physical  materials  have  properties  which  are  important  in 
relation  to  reactor  phenomena.   The  probability  that  a  material 
will  interact  with  specific  nuclear  particles  is  described 
quantitatively  in  terms  of  the  "cross  section"  for  the  reaction. 
Of  particular  interest  is  the  probability  that  an  atom  of  the 
material  will  absorb  a  neutron.   This  is  expressed  in  terms  of 

the  microscopic  absorption  cross  section  and  denoted  by  (J^  .  It 

-24   2 
carries  the  units  of  "barns";  one  "barn"  is  1C    cm  .  The  product  of 

the  number  of  atoms  in  a  unit  volume  and  the  microscopic  cross 

section  is  terrr.ed  the  macroscopic  cross  section  and  denoted  by  .2-  . 

Thus  XI*      ■  N  (J*  i  N  is  usually  expressed  in  units  of  atoms  per  cubic 
a      a 

centimeter;  therefore  2L   has  the  units  of  cm   .  The  number  of 
neutrons  in  a  given  volume  is  termed  the  neutron  density  and  denoted 
by  n.  It  is  sufficiently  accurate  to  state  that  if  the  total  neutron 
population  is  in  equilibrium,  the  neutron  density  (neutrons  per  cubic 
centimeter)  is  similarly  in  equilibrium.  Thus  <-~   ■  0.  The  partial 
derivative  is  used  since  n  is  also  a  function  of  space  coordinates. 
However,  for  a  reactor  at  or  near  critical,  the  geor  etrical  shape  of 
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the  distribution  of  n  will  not  change  significantly. 

Since  the  energy  or  velocity  of  a  neutron  is  important  in 
describing  its  behavior,  the  product  nv  where  v  is  the  neutron 
velocity  is  useful.  The  term  neutron  flux  is  defined  as  this  product 
and  is  of  necessity  a  function  of  the  energy  (or  velocity)  considered. 
This  is  denoted  by  the  symbol  $   and  carries  the  units  of  neutrons 
cm   sec   .  If  the  neutrons  are  in  thermal  equilibrium  with  the 
surrounding  medium,  they  are  called  "thermal  neutrons"  and  the  flux 
is  called  "thermal  flux".  Since  thermal  neutrons  cause  most  of  the 
fissions,  knowledge  of  cross  sections  at  thermal  energies  is  of 
considerable  importance.  Cross  section  data  are  usually  given  for 
neutrons  having  a  velocity  of  22C0  meters  per  second,  the  most 
probable  velocity  per  unit  velocity  for  thermal  neutrons  at  room 
temperature. 

If  a  reactor  is  not  critical,  there  will  be  a  net  change  in  the 
number  of  neutrons.  This  may  be  traced  to  a  combination  of  three 
causes:  leakage  out  of  the  reactor,  absorption  within  the  reactor,  and 
a  gain  froi  sources  of  neutrons.  The  leakage  per  unit  volume  is 
characterized  by  a  term  -DVp  arising  from  diffusion  theory;  D  is  the 
diffusion  constant.  The  absorption  tern;  is  written  -XI  $   utilizing 
the  definition  of  macroscopic  absorption  cross  section.  The  source 
term  may  be  written  sin. ply  as  S.   The  resulting  equation  for  the 
neutron  balance  per  unit  volume  per  second  becomes 

dvV  -  Ta  t  +  s  (1) 


at 


This  is  known  as  the  time-dependent  neutron  diffusion  equation. 
All  but  a  small  fraction  of  neutrons  froi.  fission  are  given  off 
^mediately  with  the  fissioning  of  uranium  nuclei.  These  are  termed 
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prompt  neutrons.  The  small  remaining  fraction  (C.0065  of  the  total) 
are  not  emitted  immediately  but  follow  the  radioactive  decay  of 
certain  fission  product  nuclei.  A  total  of  six  species  of  fission 

product  nuclei  which  are  precursors  for  the  emission  of  these  delayed 

235 
neutrons  are  known  for  U   .  Of  the  total  neutrons  emitted  per  unit 

tine,  the  fraction  which  is  due  to  eac:i  species  of  delayed  neutron 

precursor  is  denoted  by  fl  .,  and  A  S  2 ,  £  •   »  C.0065. 

The  remaining  possible  sources  of  neutrons  aro  from  external 

means  such  as  a  radium-beryllium  source  or  cosmic  rays. 

The  total  source  term  m-ty  be  written:   S  =  S  +  S,  +  S  for 

p    d    o 

prompt,  delayed,  and  primary  or  external  neutrons  respectively. 

The  important  characteristics  of  the  delayed  neutron  species  or 
groups  are  given  in  Table  3. 

Table  3 

Delayed  Neutron  Characteristics   1 
Group  Index  (i)    Decay  Constant,  X.  (sec   )   Fractional  Yield,  A. 

C.f 00247 

0.00138 

0.(0122 

0.00265 

C. Of  0832 

C.CCC169 


1 

0.0127 

2 

0.0317 

3 

0.115 

4 

(  .311 

5 

1.40 

6 

3.87 

Total:  ft-   0.0065C   (  -ZI^) 


It  is  important  to  note  that  the  decay  is  the  inverse  of  the 
average  lifetime  and  that  the  actual  emission  of  the  neutron  can 
occur  over  a  time  distributed  about  the  average  value.  Also  the 
fractional  yield  of  a  given  species  is  an  average  and  also  is 
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subject  to  fluctuations. 

The  discussion  of  the  time  behavior  of  a  reactor  is  greatly 
simplified  if  the  materials  in  the  reactor  core  are  homogeneous.  This 
is  approximately  the  case  with  the  AGN  2C1;  therefore  this  assumption 
is  used. 

The  technique  of  the  pile  oscillator  depends  on  the  measurer,  ent 

of  the  change  of  toe  re.-t.ctor  multiplication  factor  resulting  from  the 

oscillating  absorber  of  neutrons.  The  effective  multiplication  factor 

of  the  re.ictor,  k,  is  defined  as  the  ratio  of  the  number  of  neutrons 

present  in  the  reactor  at  a  given  time  to  the  nun  her  of  neutrons 

present  in  the  immediately  previous  generation.  The  concept  of  generation 

implies  a  characteristic  lifetime  for  the  neutrons  which,  is  the  average 

time  between  the  "birth"  and  "death"  of  a  neutron.  Since  neutrons  can 

be  lost  froi  the  reactor  either  by  absorption  in  some  nuclei  or  by 

leakage  from  the  reactor  boundaries,  it  is  useful  to  separate  the 

effects.  Thus  a  lifetime,  Jt0»  can  De  postulated  for  thermal  neutrons 

in  an  infinite  systei  from  which  there  can  be  no  leakage.  The 

reciprocal  of  the  macroscopic  absorption  cross  section  is  equal  to  the 

mean  free  path  of  the  neutron;  that  is,  the  average  distance  which  the 

neutron  trarels  until  it  is  aosorbe. .  Since  2Z  has  the  dimensions  of 

a 

cm   ,  the  mean  free  path  has  the  proper  dimensions  of  centimeters.  If 

v  is  the  characteristic  velocity  of  a  thermal  neutron,  the  neutron 

lifetime  0      is  then   — — -   .  If  a  factor  P,  which  represents  the 
Xo         rr  t 

probability  that  a  thermal  neutron  will  not  leak  from  the  reactor  is 

introduced,  the  lifetime  for  finite  geonetry  is  written  0    ■  D   P  . 

This  concept  assur.es  that  all  the  neutrons  diffusing  in  the  reactor 

are  thermal  neutrons.  A  more  correct  expression  would  include  the 
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time  a  neutron  spends  slowing  down  from  fission  energy  to  thermal 
energy.  Since  this  time  is  one  or  more  orders  of  magnitude  smaller 
than  the  thermal  lifetime,  its  neglect  will  not  be  overly  serious. 

The  factor  k  is  related  to  the  multiplication  factor  k  00    for 
an  infinite  reactor  of  the  sas- e  core  material  by  employing  an 
additional  factor  P   correcting  for  the  leakage  of  neutrons  during 
the  slowing  down  process.  Thus  k  =  k^P  P  .  The  expressions  relating 
these  factors  to  the  material  and  geonetry  of  the  reactor  are  not 
pertinent  to  this  treatment.  From  the  definition  of  k  it  is  seen  that 
if  k  =  1  the  reactor  is  critical;  that  is,  the  rate  of  creation  of 
neutrons  is  equal  to  the  rate  of  loss  of  neutrons  and  the  tiu;e  variation 
of  neutron  density  is  zero.  The  quantity  ^k   ■  k  -  1   is  a  useful 
measure  of  the  time  behavior  of  the  neutron  density.  IfAk  is  not 
zero,  the  reactor  will  not  be  in  a  steady  state  condition.  The  time 
behavior  of  the  reactor  may  also  he  related  to  the  term  "reactivity" 
defined  by  the  relation  ^  ■   ■  .  A  convenient  unit  is  the  micro-rho 
where  lMf  =  P   10   .  If  k  is  very  close  to  unity,  A  k  ^  P  . 

In  a  finite  reactor,  if  the  neutron  density  at  a  given  time  is  n, 
there  will  be  kn  neutrons  after  the  lifetime  JC   .  Thus  there  will  be  a 
net  increase  of  ZHcn ,  or  a  production  rate  of  ■  ■  *'  n  if  all  the 

neutrons  are  proi  pt.  To  correct  for  the  delayed  neutrons  the  factor 

kf3 

y.—  n,  where  ft    is  the  total  fraction  of  neutrons  which  are  delayed, 

must  be  subtracted.  Finally  there  will  be  added  those  neutrons  which 
cone  from  the  decay  of  the  delayed  neutron  precursors.  If  C.  is  the 
concentration  of  a  given  delayed  precursor  and  X   is  its 
characteristic  decay  corstant,  the 


addition  will  be  \   ^.  C 


where  the  sur.  is  over  the  six  known  delayed  neutron  groups, 
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The  e  nation  for  the  neutron  balance  is  thus 

_ .  T  .  .  _  .  +  p4  Ci  (2) 

This  equation  involves  the  unknowns  n  and  C.  A  second  equation  may 
be  written  for  the  fornation  of  the  delayed  neutron  precursors: 

dC.     .  A 

/=   kPi   n   -^.C.  (3) 

dt     — - —  1   l  y    ' 

I 
Equations  (2)  and  (3)  are  the  equations  of  the  so  called  one-group 
space-independent  model  of  reactor  kinetics. 

A  development  storting  from  the  diffusion  eouation  and  using  a 
more  sophisticated  analysis  -\12?  which  takes  some  account  of  the 
slowing  down  of  neutrons  from  fission  energies  to  thermal  energies 
n-ives  the  eouation: 


[M» -,)-»].,    ♦f^r-C*'0'    -i\ 


d(nv) 


The  factor  p  is  the  fraction  of  neutrons  which  escape  absorption 
while  slowing  down  from  fission  to  thermal  energies*  Using  the 
definitions:    k  =  k  cC   P  P„ 


J?  =4 

/-o   Vs 


o   v£; 


a 


0   -  0     P         equation  (4)  becomes 

dn   Ak      k  ft  r — ■  ,    '  ,_* 

-rr  ■  ±=L~  n r—   n   +  \    A.  C.  (5) 


i 


where   C.  =  p  PP  C . . 
i   *   f   i 

Thus  equation  (2^  will  be  in  agreeri  ent  if  the  C,  are  defined  as  the 

"effective"  concentrations  of  delayed  neutron  precursors. 
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APPENDIX  II 
DESCRIPTION  OF  THE  AGN  201  REACTOR 
t    The  AGN  201  is  a  small  low  power  reactor  having  a  licensed 
output  of  one  hundred  milliwatts.  With  increased  shielding  this  can 
be  increased  to  twenty  watts.  The  general  characteristics  of  the 
reactor  are  given  in  Table  4. 


Table  4 

Characteristics  of  the  AGN  2fl 

Power  (licensed) 

axirauni  Thermal    Flux  @    ICC   mw 
Critical    Vass 


Core  Size  (cylinder) 
Moderator  (polyethylene) 
Reflector  (graphite) 
Shield  (lead) 

(water) 
Overall  Diameter 
Overall  ileight 


ICC  milliwatts 

6 

cr 

235 


6        -2  -1 

4.5    x    1C      cm        sec 


665  gni  U  (20% 

enrichment    in   U0    ) 

10"  dia.  x  9£"  high 

11  kg 

20  cm  thick 

10  cm  thick 

55  cm  thick  (1000  gals) 

6.5  ft. 

9  ft. 


The  "glory  hole"  is  a  one-inch  aluminum  tube  which  runs 
horizontally  through  the  entire  reactor  and  passes  through  the 
center  of  the  core.  This  permits  introduction  of  various  materials 
into  the  region  of  highest  neutron  flux.  A  cross  section  view  of 
the  reactor  along  a  line  through  the  glory  hole  is  shown  in  Figure  8 

The  variation  of  thermal  neutron  flux  with  radial  position  in 
the  glory  hole  is  shown  in  Figure  9. 
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Figure  9 
Thermal  Neutron  Flux  Through  Glory  Hole 
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The  rate  of  absorption  of  thermal  neutrons  by  a  material  is 
proportional  to  the  thermal  neutron  flux  through  the  material.  The 
thermal  neutron  flux  varies  with  position  through  the  reactor.  Since 
k  is  proportional  to  the  rate  of  absorption  of  neutrons,  k  will 
depend  on  the  position  of  an  absorber  in  the  reactor.  The  relative 
effect  of  an  absorber  can  be  referred  to  its  effect  at  the  center  of 
the  reactor.  This  is  expressed  as  an  importance  function.  The 
variation  of  absorption  importance  function  with  position  in  the 
glory  hole  is  shown  in  Figure  10.  A  curve  of  squared  cosine  function 
is  fitted  to  the  importance  function  curve  over  the  diameter  of  the 
core. 
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Figure  10 
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